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Abstract. We derive gap equations for superconductivity in coexistence with 
ferromagnetism. We treat singlet states and triplet states with either equal spin 
pairing (ESP) or opposite spin pairing (OSP) states, and study the behaviour of 
these states as a function of exchange splitting. For the s-wave singlet state we find 
that our gap equations correctly reproduce the Clogston-Chandrasekhar limiting 
behaviour and the phase diagram of the Baltensperger-Sarma equation (excluding 
the FFLO region). The singlet superconducting order parameter is shown to 
be independent of exchange splitting at zero temperature, as is assumed in the 
derivation of the Clogston-Chandrasekhar limit. P-wave triplet states of the OSP 
type, behave similarly to the singlet state as a function of exchange splitting. On 
the other hand, ESP triplet states show a very different behaviour. In particular 
there is no Clogston-Chandrasekhar limiting and the superconducting critical 
temperature, To, is actually increased by exchange splitting. 
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1. Introduction 

The recent discovery of the coexistence of ferromagnetism and superconductivity in 
UGe2 P URhGe |2j and ZrZn 2 [3] has led to renewed interest in the relationship 
between ferromagnetism and superconductivity. By contrast, the relationship between 
antifcrromagnetism and superconductivity has been more thoroughly studied 4 , since 
it is relevant to many compounds, such as the cuprates borocarbides [B], heavy 
Fermion superconductors and the layered organic superconductors [S]. 

Particular interest has been focused on superconductivity on the border of a 
magnetic phase and in particular in the vicinity of a quantum critical point (QCP). 
This is observed experimentally in cuprates, several heavy Fermion systems, layered 
organics and UGe2 . It is also thought that URhGe2 may be essentially similar to UGe2 
but under the influence of 'chemical pressure' Similarly, the ferromagnetism in 
ZrZn2 also shows a QCP at high pressures. But in this case, unlike UGe2, the highest 
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superconducting transition temperatures are observed at ambient pressure, that is at 
the furthest point from the ferromagnetic-paramagnetic QCP. 

Theoretically it is thought that at or near to the QCP quantum spin fluctuations, 
can lead to spin-fluctuation induced pairing. For the case of ferromagnetic QCP this 
was first studied by Fay and Appel [5] (who also suggested that ZrZn 2 might be a 
suitable system in which to observe this effect). In this case the ferromagnetic spin 
fluctuations lead to spin-triplet pairing, by analogue with the case of superfluid 3 He. 
By contrast, in the case of quantum critical antiferromagnetic spin fluctuations spin- 
singlet d-wave pairing states are favoured J5] . 

Currently, very little is known about the superconducting pairing state in the 
ferromagnetic superconductors UGe2, URhGe and ZrZn2. If the pairing mechanism 
is indeed caused by ferromagnetic spin fluctuations, then we might expect spin-triplet 
pairing states. However, presently there is insufficient evidence in support of this 
hypothesis to be decisive. Thus it is still legitimate to consider other scenarios. 
In fact this is what we shall do here. In short, we point out that the decline of 
the superconducting transition temperature, T c , with pressure could be a simple 
consequence of p-wave pairing of arbitrary origin in an exchange field. 

In particular we will consider a simple model of the coexistence between 
ferromagnetism and superconductivity based on a parameterised electron-electron 
attractive interaction of unspecified origin. We will derive Bogoliubov-de 
Gennes (BdG) and gap equations for this model using the Hartree-Fock-Gorkov 
approximation. We will consider separately the cases of: spin singlet (s-wave) pairing, 
Opposite Spin Pairing (OSP) and Equal Spin Pairing (ESP) spin triplet (p-wave) 
states. Solving the gap equations for these pairing states, we will then illustrate some 
important properties of superconductivity in the presence of ferromagnetism. 

2. A simple model for a ferromagnetic superconductor 

We consider superconductivity arising in a Hubbard model with an effective attractive 
pairwise interaction Uij aa i , acting between electrons at crystal sites i, j with 
spins er and a'. In principle this effective interaction could arise from either 
conventional pairing mechanisms, such as electron-phonon coupling, or exchange of 
spin- fluctuations. Here we shall assume that the effective interaction is both short- 
ranged in space, namely Uij aa r ^ only for i = j or nearest neighbours, and non- 
retarded. 

In the ferromagnetic state we must also include the effective exchange field caused 
by the ferromagnetism. This enters in the model Hamiltonian as to the Zeeman 
splitting V xc . Thus the complete Hamiltonian for this model is 



where the c\J are the usual annihilation (creation) operators, hi a is the number 
operator and the o a o' are the components of the vector of Pauli matrices 



In this context we should note that the ferromagnetism of ZrZn 2 is accurately 
described by the LSDA as a weak itinerant ferromagnet. Experimentally the exchange 
splitting is clearly resolved in de Haas-van Alphen experiments ^U] an d band 




(2) 
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structure calculations (also presented in reference |1(J|) are in excellent agreement 
with these experiments. Moreover, the calculated moment (0.18/ib) is close to the 
observed moment (0.17/ub). Both the Curie temperature, Tfm, and low temperature 
magnetisation are linear functions of pressure Hence the low temperature 

magnetisation is a linear function of Tfm, in line with the predictions of the Stoner 
model. The most unusual magnetic property of ZrZn2 is that, although a field of 0.05 T 
is enough to form a single magnetic domain, the ordered moment is unsaturated up 
to 35 T El ■ This is far more naturally understood in an itinerant model such as 
LSDA or the Stoner model than, say, the Heisenberg model. On the other hand, we 
hasten to add that it is not clear whether this picture is useful for the ferromagnetic 
superconductor UGe2, since there the moments are much more strongly localised. 

Making the usual Hartree-Fock-Gorkov approximation, such that Ay CTCr < = 
—Uijo-o-i (cio-Cjtr'), and using the spin-generalised Bogoliubov-Valatin transformation, 

ha =^2u k aa'(Ri)jka' + «W (Ri )ll a , (3) 
ha' 

subject to the completeness relation 

( U k a a( R t) u kf3a(Rj) + VfeaaORiK^C-Rj)) = <%<W> ( 4 ) 

ha 

we find that the Bogoliubov de Genncs (BdG) equations for this Hamiltonian are 



/ £h + V XC 3 V xc i-iV XC 2 Att(&) 

V X ci + iV xc2 £fc - V xc3 A iT (fc) 

-Af T (-fc) -Afj(-fc) -e-k-V xc3 

\ -AJ T (-fc) -A^(-fe) -V xcl + iV xc2 

I u 1a (k) ^ 

E„{k) 



A U (fc) 

-Vxcl — iV XC 2 
-£-fe + V xc3 ) 



Ula(k) 

w To -(fc) 

V Vla(k) J 



u ia (k) 
vi ia (k) 

V Vla(k) J 



(5) 



where is the normal (that is non superconducting and non ferromagnetic) state 
energy and V xc = {V xcl , V xc2 , V xc3 ). 

The superconducting order parameter, A aa >(k), is calculated self-consistently 

from 

A aa , (k) = - i U aa , (k-q) (u aa „ (-q)v;, a „ (-g)-«: a « («)tW» (q)) (l-2/^„).(6) 

qa" 

We now introduce the Balain-Wcrthamer (BW) transformation |13H14) . 

which separates the superconducting order parameter into a singlet (scalar) part, do(k) 
and a triplet (vector) part, d(k) = (dx(k), d 2 (k), d 3 (k)). In terms of these parameters 
the BdG equations can be rewritten as 
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( £k + Vxc3 
Vxcl + iV xc 2 

-d\(k) - \d* 2 (k) 
V d* (k) + dl(k) 



-dS(fc) + ds(fc) 



-di(fe) +irf 2 (fe) 
-d Q (k)+d 3 {k) 

— £-k — V X c3 
— Vxcl + iV xc 2 



d (k)+d 3 {k) \ 

di(fc) +id 2 (fc) 

— V XC 1 — iV X c2 
— £-k + V XC 3 J 




/ M TCT (fc) \ 

u iCT (fe) 



(8) 



Using this formalism, it is also possible to calculate the free energy in the general 
case. This is given by 



At this stage one must resort either to solving these equations numerically |15H16| . 
or to studying special cases. In this paper we shall take the later approach. First, we 
begin by considering the case singlet pairing only (i.e. di(k) = ^(fe) = d 3 (k) = 0). 
In section 0] we will consider the case of only triplet pairing (i.e. when do(k) = 0). 

3. The coexistence of singlet superconductivity and ferromagnetism 

In the case of a s-wave spin singlet superconductor, it was shown by Fulde, Ferrel, 
Larkin and Ovchinnikov (FFLO) ^| ^] that the superconducting ground state 
becomes non-uniform for small external exchange fields. This solution is well known, 
and we shall not study it here. On the other hand there are also solutions which 
are spatially uniform. Whichever of these solutions is the ground state can only be 
determined by calculating the free energy for both and finding which is the lower 
solution. In strong fields the FFLO state will be the minimum, but in weaker fields 
the FFLO state will be unstable to the uniform solution. In the rest of this section 
we study the gap equations for the spatially uniform case. 

It is straightforward to show that do(k) transforms as a scalar under spin rotation. 
Thus, if there is no superconductivity in the triplet channel, we can, without loss of 
generality, rewrite the BdG equations as 




kafta 




(9) 
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(a) 




Figure 1. The four branches of the singlet spectrum in a magnetic field. Inset, 
the zero field limit where the two spin branches become degenerate. The branches 
are (a) the spectra for a =f , (b) the spectra for a =|, (c) the normal state spectra 
in zero field and (d) the singlet spectrum for V xc = 0. 



/ £fc + V xc 



V d* (k) 





-d*(k) 





E a {k) 





-d (k) 
-£-fe — V xc 


( u u (k) ^ 
u ia {k) 



d (k) 



-£- k + V xc ) 



\ ( « T «r(fc) \ 

u ia (k) 
u T<T (fe) 



(10) 



by rotating our spin reference frame so that V xc = \/V X ci 2 + V XC 2 2 + V xc ^ 2 . 

Equation (|10|l can be separated into two sets of BdG equations, so we have 



and 



d* Q (k) 
£k — Vxc 

-dS(fc) 



d (k) 

~£k + V xc 

-do(fc) 



'fc) 



wiT( fc ) 

«n( fe ) 



= £ T (fe) 



M TT (fc) 

M u( fc ) 
«u( fe ) 



(ii) 



(12) 



It is now a simple matter to regain the standard result |19| for the spectrum of a 
singlet superconductor in a spin only magnetic field: 



E a (k) 



- . LA 



\d (k)\ 2 + a\V X c 



(13) 
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with a =| = 1 and a =1 = — 1. The four corresponding energy levels are sketched in 
Fig,n Equation (|T3|) clearly reduces to the standard BCS expression for the spectrum 
of a singlet superconductor in the absence of exchange splitting as V xc — > 0. Also, 
when V xc = equations and i|12[l reduce to the usual BdG equations [20] and we 
see that we are justified in associating do(k) with the usual singlet superconducting 
order parameter A(fc). 

It is clear from (11011 that 

u CT _ CT (fe) = v aa {k) = 0, (14) 

and it can also be shown that 

u (k) = do(k) (15) 

^/(E (k)-e k ) 2 + \d (k)\2 



and 



where 



r (fc) = - E ° ik) - £k (16) 
^(E a (k)-e k ) 2 + \do(k)\* 



E (k) = V^ + Mo(fe)| 2 - (17) 

Eo(k) is, of course, of the same mathematical form as the spectrum of a singlet 
superconductor in the absence of exchange splitting. However, it is not correct to say 
that Eo(k) is the spectrum of a singlet superconductor in the absence of exchange 
splitting as the value of do(k) (although, importantly, not the value of e(k)) depends 
on V xc in general. 

Substituting our expressions for the eigenvectors of the BdG into the self- 
consistency condition © we find that the gap equation is 

d (k) = -- £ U.-Ak)^ tanh ^ 2kBT j . (18) 

In the absence of exchange splitting the gap equation regains its familiar BCS form 
|2T) . However, we note that surprisingly the exchange splitting dependence of the gap 
only enters via the Fermi (tanh) term. This means that when T = the gap equation 
becomes 

^)=-jEcW fc )|§- (19) 

ha 

which is independent of V xc . 

We must now ask what this result means physically. The most obvious conclusion 
is that, at zero temperature, the gap in independent of exchange splitting. This is 
true, but with one condition, which we will discuss below. 

The gap equation is a non-linear integral equation. And, as such, has, in general, 
more than one solution. (For example the trivial solution do(k) — is always a 
solution.) All that we have actually shown is that for any given solution do(k) is 
independent V xc at T = 0. To find the ground state we must consider all possible 
solutions and calculate the free energy of each solution. In the absence of exchange 
splitting the gap equation can be derived by minimising the free energy with respect 
to the superconducting order parameter [22] • This leads to the conclusion that the 
trivial solution is only the ground state when no other solution exists. However, no 
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Figure 2. The phase diagram of an s-wave superconductor in an exchange field 
calculated by solving the spin generalised BdG equations self consistently. Note 
that as the phase transition is first order in the presence of exchange splitting the 
free energy must be calculated for both the normal and superconducting states 
to correctly construct this phase diagram. 

such proof exists for a superconductor in a finite exchange splitting. This means that 
it is perfectly possible there to be a phase transition from the superconducting to 
normal states as the exchange splitting is increased at zero temperature. Any such 
phase transition will be 'perfectly' first order in the sense that the order parameter 
will jump from zero (above the critical exchange splitting, V^j.) to some finite value 
(below V^c) and remain at that value for all V xc < V% c . The order parameter as a 
function of exchange splitting will therefore resemble a Heaviside step function. Of 
course, as in general other superconducting phases can exist (such as the FFLO state) 
phase transitions can also occur between different superconducting phases in a similar 
manner 23. 

Such a phase transition was first studied independently by Clogston 24 and 
Chandrasekhar |2S] who both, in fact, assumed the independence of do(k) on V xc 
that we have derived above. Using this assumption they were able to show from 
simple thermodynamics that if the exchange splitting is greater than V£ c = |A(0)|/v2 
where |A(0)| is the superconducting gap at zero temperature (and zero exchange 
splitting) then the normal state has a lower energy than the s-wave superconducting 
state. This is known both of as Clogston-Chandrasekhar limiting and as Pauli- 
paramagnetic limiting. Clogston-Chandrasekhar limiting clearly applies to all singlet 
states, but does not necessarily apply to triplet states. In most superconducting 
materials hbHci < Vxc- Therefore, if a superconductor has a large upper critical 
field in comparison to the Clogston-Chandrasekhar limit this is good evidence 
for triplet superconductivity. The FFLO state can also display hbHc2 > V^.. 
Clogston-Chandrasekhar limiting has been observed in the layered organic compound 
k-(BEDT-TTF) 2 Cu(SCN) 2 [HI when a magnetic field is applied parallel to the layers 
(which prevents the formation of orbital currents due to the highly two dimensional 
nature of the material). 

To illustrate this point, we have solved the gap equation l|18(l numerically 
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for a cubic lattice. We assumed Uij aa — U J™ (i.e. an on-site interaction) 
corresponding to the case of local s-wave pairing. The comparison between the 
calculated superconducting and normal state free energies, leads to the phase diagram 
given in figure [21 This calculated phase diagram is in excellent agreement with 
that calculated from the Baltensperger-Sarma equation [23 |2H|- However, while the 
Baltensperger-Sarma equation only allows for the calculation of the superconducting- 
metal phase transition, our numerical gap equation solution allows for the evaluation 
of the order parameter at any point in T — V xc space and hence for the evaluation of 
thermodynamic variables such as the heat capacity, 

°y = E fcJpM 1 - ( E l Eka z aVxc dld ° d T )l ) ■ (20) 

and the magnetisation |19| . 



M = — t~tt V o- I d 3 fc — ^ — r-yj — ^7 . (21) 

(2tt) 3 ^ J 1 + e^te-ri/fer v ' 



where V is the volume of the first Brillouin zone. 

A numerical study of these equations |16| shows that, in an exchange field, the 
thermodynamic functions 'see' an effective gap, A e /y, i.e. 

{C v ,M, X }~e~^ (22) 

where 

A eJ/ = |A(0)| - \V XC \. (23) 



4. The coexistence of triplet superconductivity and ferromagnetism 

We will now consider the properties of a triplet superconductor in a magnetic field. 
Using a similar approach to the singlet case above we are able to derive many of the 
same physical quantities. This highlights both similarities and differences between 
the singlet and triplet cases, which may perhaps help in identifying the pairing state 
symmetry in specific ferromagnetic superconductors. 

Before we begin we will generalise a useful theorem due to de Gennes [2H. We 
begin by writing the BdG equations JSJ in a pseudo-spinor notation: 

{ -~4*_ fe -|* (fc ) ) { v_Ak) ) ~ EAk) { v_ a (k) ) ■ ^ 

Where 



£ ( , \ _ ( £k + V XC 3 V xc i — lV xc 2 \ 

i [k) ~ { V xcl + iV xc2 s k -V xc3 )> [2b) 

A - ( A TT(fc) A n (fe) \ 

=* ~ { A iT (fe) A u (fc) J' W 

^ {k) = ( Z\k) ) • (27) 
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and 

v w (k) 

Multiplying by —1, taking the complex conjugate, parity inverting and exchanging the 
rows of (|24[l leads to 



|(fc) A(fc) 

-A*(fc) -C(-k 




We have therefore shown that if 



is an eigenvector of the spin- 



as both Ea-(k) and £(fc) are even under parity inversion. 

generalised BdG equations in a magnetic field, with the corresponding eigenvalue 
( u*(—k) \ 

E a ( k ) then, — ^ ,/ I is also an eigenvector and that the corresponding eigenvalue 

is —E a (k). As a can take two values (f or J.) we have identified all of the eigenstates. 

This analysis holds for both triplet and singlet states. (For a singlet state with 
\V XC \ = it clearly reduces to the theorem of de Gennes.) The spectrum for a singlet 
superconductor in an exchange field (shown in figure ^) is clearly in agreement with 
this theorem. 

When studying triplet states, and particularly when studying the effect of 
exchange splitting on the triplet state, it is useful to introduce the notion of unitary 
and non-unitary states. For a triplet state 



4, At =i\d(k)\ 2 +i g -(d(k)xd(kr) 

and, in the absence of exchange splitting, 



(30) 



E a (k) = yjel + \d(k)\ 2 + a\d{k) x d(k)*\. (31) 

It is therefore useful to introduce the vector q(k) which is defined by 

q(k) = id(fe) X d(k)*. (32) 

It is clear that q(k) is a real vector. A unitary state is defined as any state in which 
q(k) = for all k. 

By setting the singlet order parameter, do(k), to zero we can write down the BdG 
equations for a triplet superconductor in an exchange field, 



/ 



£fc + V XC 3 

V xc \ + iV XC 2 
-dl(k) - \d* 2 {k) 
d* 3 (k) 



Vxcl — iVxc2 
£k — V XC 3 
d* 3 (k) 

d\{k) - id* 2 {k) 

( u 1a (k) ^ 
u la (k) 

V\a{k) 



E a (k) 



-di(fe) + id 2 (fe) 
da(fc) 

— £k — Vxc3 

—V xc i + iV XC 2 

( u 1a (k) ^ 
u La (k) 
w T<T (fc) 
V v l°(k) J 



d 3 (k) 
di{k)+id 2 (k) 
— Vxc\ — iV XC 2 

— £k + V XC 3 



The eigenvalues of these BdG equations are given by ^Hl 



E*(k) = \Je 2 k + ii\\V xc \ 2 + |d(fc)| 2 + o-^X(k) 



(33) 



(34) 
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A(fc) = |q(fc)| 2 +4 £ 2 Ai ||^ c | 2 + 4 A i||y, c -d(fc)| 2 +4e fc/ ^y M -<7(fc).(35) 

In zero field we clearly have the usual result Ij31|l f° r the spectrum a triplet 
superconductor . 

Again, we can also derive the expressions for thermodynamic quantities in a 
general triplet state. For example the heat capacity is given by 



Cv = Y J fka{ l BT t ] (m*? - ^\d(k)\ 2 ) (36) 
and the (vector) magnetisation, M , is given by 



M =J2 U*a(k)u la (k)f ka + V W (k)vl(k)(l - f ka ) 

k ^ 

+u\ a {k)u ]a {k)f k<J +v ia {k)v* u {k)(l - f ka ), 
-iu* u (k)u la (k)f ka - iv u (k)v* l<T (k)(l - f ka ) 
+iu| CT (fe)u T(7 (fe)/ fe<T +w ia (k)v^(k)(l - f ka ), 
+u* a {k)u ]a (k)f ka +v la (k)v* la (k)(l - f ka ) 

-u* la (k)u u (k)f ka - v l<T (k)v* a (k){l - f ka )^j . (37) 

Following the methods of Sigrist and Ueda [2H1 it can be shown ^1 that in the 
absence of exchange splitting the gap equations for a triplet superconductor are 



A al3 (k)=J2u a0 (k-k') 



4£ fcT 



d(k) 



i g(fc > x .5 fc W^ 



\Q(k)\ 



1 



g(fc) xd( fc ) 
AE ki V k ' \q{k)\ \ 2 ; 



(38) 



However, these methods do not generalise to a finite exchange splitting. Fortunately 
triplet states can be separated into three classes: those that contain only OSP states, 
those that contain only ESP states and those that contain both OSP and ESP states. 
The first two cases represent a great simplification and we will now study these special 
cases. However, it should be noted that neither of the formalisms presented below can 
deal with states that contain both OSP and ESP such as the B and B2 phases. 



4-1- Opposite spin pairing 

An OSP state is defined as any state for which d(k) x V xc = for all k. Thus, in 
this limited sense, we may describe d(k) as parallel to V xc . Much as in the case of 
singlet pairing we can, without loss of generality, rotate the system, recalling the d(k) 
transforms as a vector under rotation, so that 
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/ £k + V xc 





£k - V xc 

dm 





EJk) 



(fe(fc) \ / u u (k) \ 

£-k-V xc WTcr(fc) 

-e-h + T4c / V v i°( k ) J 

/ Ut<r(fe) \ 
•fe) 
fe) 

V v ia {k) J 

Again we can separate (|39[1 into two BdG equations and hence, in a similar manner 
to which we derived the singlet gap equation, we find that the gap equations for OSP 
triplet superconductivity are 



(39) 



*(*) = - 1 - V U„ (k)^\ tanh (BoW+'V™) 



(40) 



Note that this equation is of precisely the same mathematical form as the singlet gap 
equation (|18fl . Both the phase diagram and the effective gap 'seen' by thermodynamic 
probes are the same as we earlier found for singlet superconductivity. However, this 
time the effective gap 'seen' by thermodynamic probes is given by ^3 El 

(41) 



A eff = \d(k F )\ - \ V X , 



where |d(fej?)| is the mean gap at the Fermi surface. 

All singlet pairing states are, by definition, OSP states. Thus it appears that, in 
the presence of exchange splitting, the important property of a state is whether it is 
an OSP or an ESP state, not whether it is a triplet or a singlet state. 

4-2. Equal spin pairing 

An ESP state is defined as any state for which d(k) ■ V xc = for all fe. Thus, in 
this limited sense, we may describe d(k) as perpendicular to ~V XC . In this case, for 
V xc = (0, 0, — V xc ), the spin triplet BdG equations are 



/ £fc — V xc 


-A| T (-fe) 

V o 





£fc + V X c 


-Ah(-fc) 
= E a {k) 



A TT (fc) 


-£_fc + V xc 







A U (fe) 


-£-k — V xc ) 



\ t «T«r(*) \ 

u la (k) 
v u (k) 

V «l«r(*) / 



(42) 



/ u w (k) \ 
u ia (k) 

V]a{k) 

\ v ia {k) J 

We can now easily separate the BdG equations into a pair of BdG equations for 
up electrons, 



e k - V xc A TT (fe) 
-A| T (-fc) -e-fc + 14, 

and a set of BdG equations for down electrons, 



u To -(fc) 
v 1a (k) 



= E a (k) 



U<[a{k) 

v^(k) 



(43) 
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- k 4/-k -tf-l )(:;:$ )-*«(.,.<»)<«> 

Using the self-consistency condition © we easily find that the gap equations are 

A ^ k) = _ ^ u aa {k-k>)£ {k > ) {i 2f ^ y (45) 

k' 

with 



Eka = ]f (e k - aV xc ) 2 + \A aa (k)\i (46) 

As T — > T c from below, |4 | -> and hence J5 ff (fe) -> e(fc) + aV xc . Therefore 
the gap equation becomes 

ffffV ' ^2(e(fc')-a^ c ) V 2fc B T ^ ff<A ^ V j 

Thus, near Tc the gap equation is linear. This allows Tc to be determined very 
accurately. Further by comparing the transition temperatures of various symmetries 
one can find which has the highest transition temperature and hence which state 
occurs for T<T C . 

Clearly, one cannot, in general, use the linearised gap equation to study transitions 
from one superconducting state to another as the gap equation can no longer be 
linearised below the first superconducting transition. The exception to this rule is the 
transition from an ESP state with only one type of pairing to an ESP state with both 
j'f and H pairing (an example of such a transition is the transition from the Ai phase 
to the A2 phase), because of the complete separation of the spin- up and spin-down 
subsystems in the presence of exchange splitting and the absence of opposite spin 
pairing or spin flip processes. 

We solved the linearised gap equations l|47|l numerically for parameters chosen 
of ZrZn 2 (see j3U] for a discussion). To do this we used a simple cubic tight binding 
model and a k-space integration mesh of 10 9 points. We use such a large array for 
two reasons. A fine integration mesh is required to accurately determine the density 
of states (DOS). Our method (implicitly) requires an accurate calculation of the spin 
dependant DOS, D a {ep)- This is particularly important in our case as we are varying 
the exchange splitting and thus we are changing the D ty {eF) 1 so any errors in evaluating 
D a (e) will lead to significant errors in our calculation of the variation of Tc with V xc . 

We show the results of our numerical calculations in figure |31 The line is a cubic 
curve fitted to the numerical data. For any given exchange splitting, V xc , there are 
two transition temperatures, corresponding to the two separate spin components of 
the ESP order parameter. We have plotted the transition temperature for ft pairing 
on the positive V xc side of the graph and the transition temperature for J, J, paring on 
the negative V xc scale. There are several reasons for plotting the data in this way. 

(i) In this way the graph shows the behaviour of the ft pairing state over a full 
range of exchange splitting, from positive to negative. 

(ii) We see that the point V xc — is not a special case, and the curve is smooth 




Figure 3. The results of our numerical solution of the linearised gap equations 
are shown by the points. The line is a fit to the calculated points by a cubic 
equation. 



there. 

(iii) We also have a larger data range to fit over, and thus increase the accuracy 
of the cubic fit. 

Zero exchange splitting is not a special point because in both the non-linear and 
linearised gap equations exchange splitting is mathematically equivalent to a change 
in chemical potential. Thus, the graph plotted in the manner shown in Fig. [3] can 
also be interpreted as a plot of critical temperature of the of the triplet A phase as a 
function of the chemical potential in zero exchange splitting. 

We now plot the critical temperature for both and || pairing on the same 
graph (figure HJ. This plot shows is then the (V xc , T) superconducting phase diagram 
for our model. (This, of course, assumes that no further phase transitions occur at low 
temperatures.) The higher transition temperature is the transition to the Ai phase 
(where only ft pairing occurs) and the second transition is a transition to the A2 
phase (where || pairing begins). In the paramagnetic state (the line V xc = 0) the 
superconducting state is an A phase as the superconducting order parameter is the 
same for both the ft an d II pairing states. (The A 2 phase becomes the A phase via 
a cross over, rather than a phase transition.) 

The phase diagram shown in figure0]is clearly equivalent to the A1-A2 splitting of 
3 He in a magnetic field. Experimental measurement of this phase transition in 3 He due 
to Remeijer et al are reported reference At first sight figure 0] and reference jSJ 
appear rather different, however the are in fact almost identical, as we will now show. 
The dimensionless measure of the exchange splitting for the Remijer et al experiments 
is , where Tp is the Fermi temperature and /i„ is the nuclear magnetron for 3 He, 
while for our calculation the dimensionless exchange splitting is given by ^jjf- where 
W — 16t is the bandwidth. The experiments of Remeijer et al were not performed at 



Th 




constant pressure, which complicates the analysis somewhat, however they conclude 
that 



TA \k B T F J ' " \k B T F 



(48) 



where a = 36.3 ± 0.91 and b = 522 ± 17 in the range < p^- < 0.01 at an 

° — k B If — 

effective pressure of 3.4 MPa i.e the splitting is, to a very good approximation linear. 
The equivalent exchange splitting in our calculations is V xc < 0.01 W = 0.01 eV. It 
can clearly be seen from figure 0] that our calculations give an approximately linear 
splitting between the Ai and A2 phase transitions over the range of exchange splitting 
< V X c < 0.01 eV. Hence our results are consistent with the what is known about 3 He. 
(Although, of course, we had no right to expect this agreement as our parameters where 
chosen for ZrZn2 and not 3 He.) Further this illustrates the fact that ferromagnetic 
superconductors will provide an excellent laboratory in which to study the splitting of 
the Ai and A2 phase transitions (and the non-linear splitting in particular) over a far 
greater range of exchange splitting than is possible in 3 He. Further, when the effects 
of scattering from non-magnetic impurities are included this model gives results that 
are qualitatively consistent with the observed pressure dependence of Tq in ZrZn2 

5. Discussion 



We have derived gap equations for superconductivity in coexistence with 
ferromagnetism. We have done this for s-wave singlet states and for p-wave triplet 
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states with either ESP or OSP pairing. We used these gap equations to study the 
behaviour of these states as a function of exchange splitting. 

For the singlet state we found that our gap equations reproduced the Clogston- 
Chandrasekhar limiting behaviour and the phase diagram of the Baltensperger-Sarma 
equation (neglecting the possibility of an FFLO state). We also showed that the 
singlet gap equation leads to the result that the superconducting order parameter is 
independent of exchange splitting at zero temperature. This fact was assumed in the 
derivation of the Clogston-Chandrasekhar limit. 

OSP triplet states showed a very similar behaviour to the singlet state in the 
presence of exchange splitting. This leads to the conclusion that the effect of exchange 
splitting on a superconducting state is determined by whether the state contains OSP 
or ESP. (All singlet states are, by definition, OSP states.) 

In contrast, ESP triplet states show a very different behaviour in an exchange 
field. In particular there is no Clogston-Chandrasekhar limiting. Further, Tq is 
actually increased by exchange splitting because D a (eF) is changed by the exchange 
splitting and Tc is dependent on D CT (£j?). This effect is well known in 3 He, but has 
previously only been studied in a Ginzburg-Landau formalism [22] ■ The gap equations 
presented here will allow for far more detailed study of both the increase of Tq< and 
for the study of the splitting of the Ai and A 2 phases by exchange splitting. 

If the experimentally occurring ferromagnetic superconductors are ESP triplet 
pairing states, as seems likely from the absence of Clogston-Chandrasekhar limiting, 
then these systems will allow for study of this effect at far greater exchange splittings 
than can be archived with magnetic fields in 3 He. The gap equations presented here 
will also be useful for studying these materials in their own right, in particular we 
hope that the will prove useful for identifying the superconducting pairing symmetry 
of these ferromagnetic superconductors. Our formalism is quite general, and can be 
applied to more realistic band structures and pairing models, although the additional 
complication of the vector potential will have to be overcome before one can make 
complete theoretical predictions for these materials. 
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